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ABSTRACT 

There is a striking and unexplained dearth of brown dwarf companions in close orbits (< 3AU) around stars 
more massive than the Sun, in stark contrast to the frequency of stellar and planetary companions. Although 
rare and relatively short-lived, these systems leave detectable evolutionary end points in the form of white 
dwarf - brown dwarf binaries and these remnants can offer unique insights into the births and deaths of their 
parent systems. We present the discovery of a close (orbital separation ~ 0.006 AU) substellar companion to 
a massive white dwarf member of the Praesepe star cluster. Using the cluster age and the mass of the white 
dwarf we constrain the mass of the white dwarf progenitor star to lie in the range 3.5 - 3.7 M Q (B9). The high 
mass of the white dwarf means the substellar companion must have been engulfed by the B star's envelope 
while it was on the late asymptotic giant branch (AGB). Hence, the initial separation of the system was ~ 2 
AU, with common envelope evolution reducing the separation to its current value. The initial and final orbital 
separations allow us to constrain the combination of the common envelope efficiency (a) and binding energy 
parameters (A) for the AGB star to a\ ~ 3. We examine the various formation scenarios and conclude that the 
substellar object was most likely to have been captured by the white dwarf progenitor early in the life of the 
cluster, rather than forming in situ. 

Subject headings: binaries: close — brown dwarfs — white dwarfs 



1. INTRODUCTION 

There is a known dearth of brown dwarf companions to 
solar-type stars with orbital periods <5 years (equivalent to 
orbital separations <3AU) when compared with lower mass 
planetary companions or mo re massive stellar companions 
dGrether & Lineweav er 2006). There is also some evidence 
that th is paucity of objects may exte nd to much larger separa- 
tions (McCarthy & Zuckerman 2004). The reason for the lack 
of brown dwarf companions at these separations is unknown, 
but it is likely related to the formation mechanisms involved. 

The difficulties in identifying brown dwarfs with early type- 
companions mean the most extreme examples of these bi- 
nary systems are found in a highly evolved form: white 
dwarf - brown dwarf binaries. However, detached brown 
dwarf and very low-mass stellar companions to white dwarfs 
are rare; the fractio n of L-type secon daries to white dwarfs 
is just 0.4±0.3% dSteele et alj 120111) . Proper motion sur- 
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veys and searches for IR exces ses have so far found only a 
handful of confirmed examp les (Becklin & Zuckerman ll 1 998k 
Farihi & Christ opher! 12004 iMaxted et al.ll2006t ISteele et all 
Burleighet all 1201 U iDav- Jones et alj l201lT 
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none of which have a reliably determined 
age independent of the white dwarf p arameters. Only in 2 sys- 
tems, WD0137-349B (L8, 0.053 M^:lMaxted et alJl200 6) and 
GD1400B (L6-L7, 07-0.08M w ; Farihi & ChristopheEool 
iBurleigh et al.ll201lb (P or b = 116 minutes and 9.8 hours re- 
spectively) is the brown dwarf known to have survived a phase 
of common envelope (CE) evolution. This phase of binary star 
evolution involves the brown dwarf being engulfed by, and 
immersed in, the expanding atmosphere of the white dwarf 
pro genitor as it evolves away fro m the main sequence (see 
e.g. iDavis. Kolb & K nigge 20121). In this letter we present 
the discovery a new white dwarf - substellar binary system in 
the Prasespe open star cluster. We show how the cluster age 
along with the mass and cooling age of the white dwarf can be 
used to place two independent limits on the mass of the white 
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FIG. 1 . — Radial velocity phase diagram: The radial velocity data folded on 
the most likely period of 4.2 hours. v*sin i= 11.31±1.55 kms -1 . 

dwarf progenitor star and, additionally, to constrain the initial 
orbital radius of the brown dwarf. We use this information 
to examine the physics of common envelope evolution and to 
test formation models for the original system. 

2. OBSERVATIONS 

We originally obtained high resolution optical spectra of 
the 0.798±0.006 M white dwarf WD0837+185 to confirm 
its membership of the 625 ±50 Myr old Praesepe open star 
cluster and discovere d that its radial velocity was varying 
dCasewell et al.1 120091) . but that it was not a double lined 
spectroscopic binary. Subsequently, 22 follow-up observa- 
tions were acquired between 2008/02/07 and 2008/03/1 1 from 
the Very Large Telescope's Ultraviolet Echelle Spectrograph. 
These data were obtained with exposure times of 20 min in 
thin cirrus, with seeing between 0.5 and 2.0". Each pair of 
consecutive datasets (obtained ~ 1 minute apart) were coad- 
ded to increase the S/N. The data were obtaine d with the same 
grating settings, reduced and analysed as in ICasewell et al.l 
(2009). These new measurements verified the variation with a 
best fitting period of 4.2 hours, confirming the presence of an 
unseen companion. The velocity semi-amplitude calculated 
from the system parameters is, K = v* sin;= 1 1 .3 1 ±1 .55 
km s" 1 (Figure [TJ, giving a minimum mass for the compan- 
ion M sin/ w 25 Mj up . We also obtained photometry 
from the United Kingdom Infrared Telescope (UKIRT) us- 
ing the UKIRT Fast Track Imager in the J, H, and K, bands. 
These data had exposure times of 3000s in the K, 1200s in 
the H and 600s in the J band using a 5 point jitter pattern. 
The data were reduced using the STARLINK package ORAC- 
DR, and the photometry performed using IRAF. We also ob- 
served WD0837+185 with the Spitzer space telescope us- 
ing IRAC in the [3.5] and [4.5] micron bands (Cycle 7, Pro- 
gramme ID:57771, PI: Casewell). These data were reduced 
using the MOPEX pipeline, and the aperture photometry per- 
formed using APEX before the pixel phase, array location and 
aperture corrections were applied. 

3. RESULTS 

A comparison of optical photometry from the Sloan Dig- 
ital Sky Survey (SDSS), near infrared photometry obtained 
with UKIRT, the UKIRT Infrared Deep Sky Survey, and mid- 
infrared images from Spitzer with a pure Hydrogen atmo- 
sphere whit e dwarf model for WD 0837+185 (T eff = 15000 K, 
log g= 8.3; ICasewell et al.l 120091) . and the Hydrogen model 
combined with observed spectra of a T5, and a T8 dwarf 
showed that no excess emission, indicative of a companion, 
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FIG. 2.— Photometry and models of WD0837+185: SDSS ugriz (trian- 
gles), UKIRT ZYJHK (boxes) and Spitzer [3.6], [4.5] (diamonds) magni- 
tudes shown with a DA white dwarf model spectrum (T c ff= 15000 K, log 
g=8.3). 3cr errorbars are shown on the Spitzer datapoints. A DA+T5 com- 
posite spectrum is also shown as a dashed red line, and a DA+T8 spectrum 
as the dotted blue line. The T dwarf spectra are rea l data and are the objects 
2MASSJ 05591914-14044 88 iCushing et al., 2006) and 2MASSJ0415I954- 
0935066 I Saumon et al., 2007), but the spectra are not continuous as they are 
M and L band spectra. There are gaps between 4.1 and 5.2 fim. See the 
electronic edition of the Journal for a color version of this figure. 

is seen. However, if we increase the errors to the 3<r level 
there is a possibility of an excess in the [3.6] and [4.5] micron 
wavebands if the white dwarf has a T8 or cooler companion. 
(Figure[2]i. WD0837+185 has a luminosity, proper motion, ra- 
dial velocity, cooling age and gravit ational redshift con sistent 
with being a member of Praesepe (Casewel l et al.ll2009h . If 
the unseen companion were another white dwarf, this would 
increase the total luminosity and make it inconsistent with 
Praesepe membership. Even a high mass (~ 1.38M Q ) and 
therefore small radius white dwarf with a cooling age ~ the 
cluster age (giving an effective temperat ure of 25000K) would 
be detected in the optical photometry dHolberg & Bergeron 
l2006t iTremblav. Bergeron & Gianninasll201 ll) . Moreover, if 
the secondary were another white dwarf or even a neutron 
star the radial velocity solution would require an extremely 
low inclination orbit (^1-4 degrees; Probability < 2.4xl0~ 03 ) 
to hide its influence. Hence, we rule out the possibility 
that WD0837+185B is another white dwarf or a neutron star. 
WD0837+185 is not coincident with an X-ray source in the 
ROSAT all sky survey, eliminating the possibility the system 
contains an accreting black hole. The probability of the sys- 
tem containing a non-accreting black hole is 6xl0~ 6 , as the 
inclination would need to be extremely low (~0.2 degrees) to 
match radial velocity measurements. 

The photometry in the [4.5] micron band gives the maxi- 
mum possible mass of a non-degenerate secondary as 30Mj up 
dBaraffe et al.l2003l) . Combined with the minimum mass from 
the radial velocity solution (M sin i « 25 Mj up ) this strongly 
suggests that the companion is substellar and is likely to be 
a methane atmosphere T-type brown dwa rf (T8 or later) with 
an effective temperature of «900- 1 1 00 K dBaraffe et alj2003l) 
(Figure |2j. 

We have also obtained V band photometry obtained ev- 
ery 30s over a total of 5 hours from the Isaac Newton Tele- 
scope on La Palma in March 2009 to investigate the possi- 
bility of an eclipse in the system. The data show variation 
at the 0.74% level (peak to peak) on the orbital timescale 
(Figure |3J. This may be indicative of irradiation by the 
high ultraviolet flux of the white dwarf, possibly leading to 
substantial temperature differences between the "day" and 
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latter being defined by 



FIG. 3. — Light curve from March 2009 binned by a factor of 7 folded on 
the orbital period of 4.2 hours. The peak to peak variation is at 0.74%. 

"night" hemisphe res of the T dwarf atmosphere as is seen 
in WD0137-349 dMaxted et al.1 120061) . No eclipse is seen, 
but the probability of eclipse in this syste m is only ^9% an d 
any eclipse would only last ^3.5 minutes dFaedi et alj|201 lb . 
which is difficult to detect with our sampling frequency (30s 
exposure + 30s readout). The main-sequence lifetime of the 
white dwarf progenitor is constrained by the difference be- 
tween the cooling age of WD0837+185 (313±5 Myr, deter- 
mined from the white dwarf temperature, gravity and appro- 
priate cooling model; Fontaine, Brassard, & Ber geronll2001l) 
and the age of the Praesepe cluster (625±50Myr), which 
limits the progenitor mass to 3.48±O.23M equivalent to a 
spectral type B9. Intriguingly, the orbit of WD0837+185B 
(orbital separation —0.006 AU ~1.24R ) is w ell within the 
main sequence radius of a B9 star (R ~ 3R ; Mari go et alj 
120081) . WD0837+185B is unlikely to have been captured 
into this current orbit as the star-crossing timescale in Prae- 
sepe is a few Myr and the low stellar density (20 starspc" 2 ; 
Kraus & Hillenbrand 2007) makes the typical close encounter 
(<SC 1AU) time-scale much longer than the age of the Uni- 
verse. The system must therefore have gone through a phase 
of CE evolution. This brief evolutionary phase is a key stage 
in the formation of short period binary systems and occurs 
when the lower mass companion becomes engulfed by the 
nuclear-evolution driven expansion of the primary star's en- 
velope. The binary orbit then decays rapidly as drag forces 
unbind the primary's envelope at the expense of orbital en- 
ergy. 

The mass of WD0837+185 (~O.8M ) is at the very up- 
per end of core masses attainable within the errors on the 
progenitor star mass. Models relating initial stellar mass to 
core mass (taken to be the white dwarf m ass) at first ther- 
mal pulse (Karakas, Lattanzio & Pols 2002) give a lower limit 
of ~ 3.6M Q for the progenitor mass, while the cooling age- 
cluster age argument above provides an upper limit of ~ 
3.7M Q . The very limited overlap between these indepen- 
dent mass estimates show that the progenitor star engulfed 
the brown dwarf very late in its evolution on the AGB. The 
radius of the progenitor at this point, and hence the orbital 
separation of WD0837+ 185B at the onset of the CE phase, 
would have been ~ 2AU dVentura & Marigo 2009). This esti- 
mate places the orbital separation of the original system well 
within the region where there is an observed dearth of brown 
dwarf companions to solar-type stars. 

Treatments of CE evolution are usually parameterized in 
terms of the efficiency with which orbital energy is trans- 
ferred to the envelope of the primary star (a) and a param- 
eter governing the binding energy of the primary's envelope 



^bind — " 



(1) 



where £bind is the binding energy, Mi is the progenitor (pri- 
mary) mass, Menv is progenitor envelope mass and Ri is the 
progenitor's radius. Many population synthesis calculations 
treat A as a constant, but its value is poorly constrained. Treat- 
ments differ in the inclusion (or not) of the internal energy 
of the stellar matter. We are able to place explicit limits 
on A because we know that the CE phase began when the 
WD0837+185 progenitor was on the late AGB, which fixes 
the original orbital radius as well as the core and the envelope 
mass of the progenitor. Knowledge of the initial (a, — Ri) and 
final (a/) orbital separations allows the combination aX to be 
extracted from the energy balance equation 



aX 



GM core Mi GM X M 2 \ GMiM, 
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where M2 is the brown dwarf mass and M core (= Mi - M env ) 
is the progenitor core mass. For a progenitor mass 3.6 M 
and a companion mass of 30Mj up the combination aX —3 is 
required to place the brown dwarf in its current orbit from 
an initial orbit — 2AU, assuming no significant mass loss by 
the progenitor at the point of contact, (any mass loss would 
cause a proportionally lower estimate of aA). In the case of 
maximum efficiency (a =1), A —3 is very low for a highly 
evolved 3.6 M star if the internal energy of the stellar mate- 
rial is taken into account when calculating A, but is in reason- 
able agreement with calculation s including only gravitational 
binding energy dXu&Lill2010l) . Highly evolved, late AGB 
primaries of this mass are precisely those predicted to have 
the highest values of A when including the contribution of the 
internal energy of the stellar matter, and a low A for these ob- 
jects is a strong hint that this parameter is of order unity for 
all primaries. 

4. INTERPRETATION 

WD0837+185B may have formed (at its original orbital 
separation) in one of two ways: in a manner similar to So- 
lar System and extra-solar giant planets, or as an extreme 
mass-ratio binary star. The planetary channel assumes for- 
mation in a disk around the newly-formed star, and con- 
sists of two options: core accretion or gravitational instabil- 
ity. Core accretion is not a promising formation mechanism 
here: tidal torques from the planet strongly suppress gas ac- 
cretion from the disk on to the p lanet for masses > 5Mj up 
dD'Angelo. Henning & Klevll2002l) . and it is unlikely that ob- 
jects as massive as as ~20Mj up can ever form in this manner. 
Moreover, the growth time-scale for such a massive object is 
at least comparable to the typical lifetime of protoplanetary 
disks around B stars (w hich are estimated to be Myr or less; 
Hillenbrand etlfl[l992h . 

Disk fragmentation via gravitational instability can re- 
sult in the formation of much more massive objects, and 
is a plausible scenario for brown dwa rf formation (e.g., 
Stamatellos, Hubber & Whitworthl 120071) . but carries the 
caveat that protostellar disks are gravitationally unstable only 
at large radii. For a 3.5 M star a self-luminous disk is only 
unsta ble at radii >70 AU dMatzner & Levinl 120051; iRafikovl 
2009), and this critical radius increases when irradiation from 
the star is taken into account. The question then becomes 
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whether a brown dwarf can migrate inwards to 1 -2 AU from its 
formation radius at ^100 AU. Dynamical "hardening" of the 
system via repeated encounters with other cluster stars is ex- 
tremely unlikely, as repeated such interactions are more likely 
to disrupt the system than shrink its orbit; also, the low stel- 
lar density in Praesepe essentially rules out this mechanism. 
Inward migration via gravitational interactions with other 
brown dwarf- or planetary-mass objects (so-called planet- 
planet scattering) is possible, but again unlikely: the process 
is chaotic, but simulations find that inward migration via this 
mechanism is usually modest (e.g., Raymond et alJ2008l) : the 
probability of being scattered from ^100 AU to ~l-2 AU is 
very small. The final possibility is gas-driven migration be- 
fore the dispersal of the protostellar disk. Giant planets typ- 
ically migrate in the Type II regime. Although the discov- 
ery of two hot Jupi ters around stellar members of Praesepe by 
lOuinn et"afl (120121) shows that planet formation and migration 
in the Type II regime has occurred in this open cluster, it is 
unlikely that objects as massive as ^25 Mj up can ever form in 
this manner. For an object of 25-30 Mj up Type II migration is 
strongly suppressed and proceed s on a time-scale muc h longer 
than the lifetime of the gas disk (ISyer & Clarkdl 19951) . but re- 
cent simulations have found that migratio n in gravitationally 
unstable disks can in fact be very rapid (|Cha & Nav akshin 
1201 U iBaruteau. Meru & Paardekooperil201 ll) . It is unclear, 
however, whether it is possible to halt the rapid migration of 
20-30 Mj up objects at ~1 AU, and prevent them falling all the 
way on to the central star. This mechanism cannot be ruled 
out without further investigation, but we do not consider it to 
be the most likely formation channel for WD0837+185B. 

Binary stars form from the fragmentation of star-forming 
molecular cores, but are affected by the cluster dynamics as 
they evolve. Extreme mass-ratio binaries in close orbits such 
as this one (q ~ 0.01) are very rare. Forming such systems 
in a similar manner to binary stars is challenging, because 
circumstellar material preferentially accretes on to the sec- 
ondary, increas ing its mass and dri ving the binary mass-ratio 
towards q = 1 ([Artymowicz 1983). Many theoretical stud- 
ies have investigated how dynamical interactions may harden 
initially wide binaries, or destroy them, but the fundamental 
phy sics of binary formati on remains poorly understood (see, 
e.g.. lGoodwin et afll2007l and references therein). 

The one remaining formation theory left for consideration 
is thus dyn amical capture. Recent numerical simulations 
dBatdl201ll) show that in a N~ 100-500 cluster a few objects 
are formed this way, though due to the low number statistics 
it is not clear how frequently extreme mass-ratio objects are 
captured into ~ AU orbits. Praesepe has more membe rs than 
in this example (N^IOOO; iKraus & Hillenbrand! 120071) . and 
at its current age, is likely to have undergone dynamical evo- 
lution, decreasing the cluster population by as many as half, 
and preferentially ejecting the lowest mass cluster m embers 
(Ide La Fuente Marco s & de La Fuen te Marcos] 120021) . The 



most plausible formation scenario for WD0837+185B is 
therefore likely to be dynamical capture where the brown 
dwarf has been inserted into the binary after its formation. 
This is likely to be a common formation scenario for those 
relatively rare oases in the brown dwarf desert: high mass ra- 
tio main sequence plus brown dwarf pairs found b y radial ve- 
locity surveys with separations of a few AU (e.g.. lOmiva et alj 
1201 



5. CONCLUSIONS 

We confirm that WD0837+185 is a radial velocity vari- 
able object and conclude from optical, near- and mid-IR pho- 
tometry that the probable companion is a 25-30Mj up late T 
dwarf. Optical photometry also tentatively indiactes that the 
white dwarf is irradiating its substellar companion, although 
no eclipse is seen in the data. 

Using the cluster age and the mass of the white dwarf we 
constrain the mass of the white dwarf progenitor star to lie in 
the range 3.5 - 3.7 M Q (B9). The high mass of the white dwarf 
means the substellar companion must have been engulfed by 
the B star's envelope while it was on the late AGB. Hence, the 
initial separation of the system was ^2 AU, with common en- 
velope evolution reducing the separation to its current value. 
The initial and final orbital separations allow us to constrain 
the combination of the common envelope efficiency (a) and 
binding energy parameters (A) for the AGB star to aX ~3. We 
examine the various formation scenarios and conclude that the 
substellar object was most likely to have been captured by the 
white dwarf progenitor early in the life of the cluster, rather 
than forming in situ. 
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